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ABSTRACT 

We study the effect of quasar feedback on distributions of baryons and properties of 
intracluster medium in galaxy groups using high-resolution numerical simulations. We use 
the entropy-conserving Gadget code that includes gas cooling and star formation, modified 
to include a physically-based model of quasar feedback. For a sample of ten galaxy group- 
sized dark matter halos with masses in the range of 1 to 5 x 10^^ Mq / h, star formation is 
suppressed by more than 50% in the inner regions due to the additional pressure support by 
quasar feedback, while gas is driven from the inner region towards the outer region of the 
halos. As a result, the average gas density is 50% lower in the inner region and 10% higher in 
the outer region in the simulation, compared to a similar simulation with no quasar feedback. 
Gas pressure is lowered by about 40% in the inner region and higher in the outer region, 
while temperature and entropy are enhanced in the inner region by about 20-40%. The total 
group gas fraction in the two simulations generally differs by less than 10%. We also find a 
small change of the total thermal Sunyaev-Zeldovich distortion, leading to 10% changes in 
the microwave angular power spectrum at angular scales below two arcminutes. 
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1 INTRODUCTION 

Galaxy clusters and groups are the largest gravitationally bound 
objects in the universe, and they dominate the total baryon con- 
tent of the universe. Their spatial distribution and mass function 
contain information about the formation and evolution of large- 
scale structure, which in turn constrain a variety of fundamen- 
tal cosmological properties including normalization of the mat- 
ter power spectrum, the cosmic baryon density, and dark mat- 
ter properties. However, in order to use them as a cosmological 
probe, it is necessary to understand their astrophysical properties, 
and in particular their baryon physics. This issue is of particular 
current interest due to upcoming arcminute-resolution microwave 
sky survey s like ACT I Kosowskv et al. 20061 : iFowler et al. 2007h 
and SPT jRuhl et al 2004h . which will image galaxy clusters 
via the Sunyaev-Zeldovich distortions to the cosmic microwave 
blackbody spectrum from the hot electrons in the cluster gas 
dSunvaev & Zeldovich 1980h . 
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The majority of baryons in clusters and groups are in the 
form of hot intracluster gas rather than than individual galax- 
ies. Properties of the Intracluster Medium (ICM) have been 
studied through a c ombination of X-r a y and radio observ ations 
jNulsen et al. 20051 : iHeinz et al. 20021 : iFabian et al. 2000l) . Al- 
though the dark matter distribution in galaxy clusters follow a self - 



similar relation jPointecouteau et al. 20051: Jvikhlinin et al. 2006 ). 
the hot gas does not jSanderson et al. 20031 : IPopesso et al. 2005 ). 
Additional non-gravitational sources of heating are required to ex- 
plain the observations. One interesting and plausible possibility is 
the energy radiated from quas ars or Active Galactic Nuclei (AGN) 
and deposited into the ICM jKaiser 199ll: Ivalageas & Silk 199S : 



Nath & Rovchowdhurv 20021 : 



Thacker et al. 20061) . whichwe study in this work. 



Scannapieco et al. 2005 



The best arena in which to study the impact of various feed- 
back mechanisms is galaxy groups. Massive clusters with deeper 
gravitational potential wells are likely to have their global thermo- 
dynamic and morphological properties less affected by feedback. 
In comparison, galaxy groups have shallower potential wells while 
still having enough gas to display the effect of feedback on the 
ICM. Galaxy groups have recently been obs erved in X-rays at red- 
shifts as large as z = 0.6 jWillis et al. 2005 ). In the optical band, 
jTago et al. 2007t) have compiled group catalogs from the SDSS 
Data Release 5 catalog. Evidence for heating by a central AGN 
or radio source in galaxy groups and clusters has been the subject 
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of several recent papers ([Croston. Hardcastle & Birkinshaw 20051 : 
Ijetha et al. 20071 : 1 Sanderson et al. 20051) . These observations show 
excess entropy in cluster cores, which suggests that some heating 
process must act to offset cooling. 

In recent years, cosmological simulations including dark 
matter and gas have been able to follow the evolution of 
individual galaxy groups and clusters. A number of studies 
have investigated the cluster baryon fraction and its evolu- 
tion in numerical simulations. Adiabatic simulations that do 
not include radiative cooling find cluster bar yon fractions 
around 0.85 of the universal baryon fracti on tevrard 19901; 



Metzler & Evrard 1994 



Navarro. Fre nk & White 1995 





Lubin et al. 1996; Eke, Navarro & Frenk 1998; Frenk et al. 1999 
Mohr, Mathiesen & Evrard 1999 ; Bialek. Evrard & Mohr 2001) . 
Preheating the gas r educes the fraction further 
teialek. Evrar d & Mohr iooH: [Borgani e t al. 20021 : 



[Muanwong et al. 20021 : Kav. Thomas & Theu ns 2003). When 



cooling, star formation and other feedback processes are included, 
the baryon fraction is higher than th at obtained from adiabatic sim- 
ulations I Muanwong et al. 20021: Kav. Thomas Sl Theuns 20031 : 
IValdamini 20031 : lEttori et al. 20041 : iNagai et al. 2007h . This leads 
to an "overcooling" problem and indicates an additional feedback 
mechanism. 

In the current study, we analyze the effect of quasar feedback 
on the baryon distribution and thermodynamics of hot gas in galaxy 
groups. We also study its implication for the Sunyaev-Zeldovich 
angular power spectrum, which receives a dominant con tribution 
from high-redshift halos. Ref. jKomatsu & Seliak 2002h showed 
that the thermal SZ angular power spectrum provides a strong 
constraint on the normalization of the matter power spectrum, 
as. Upcoming SZ surveys like ACT or SPT will have sufficient 
sensitivity to determine with an accuracy limited by uncertainty 
in the theoretical model. Also, the kinematic SZ effect is a 
measure of bulk motions in the universe and may be a co mpetitive 
probe for studying cosm ology (Sehgal, kosowskv & Ho lder 20051: 



Bhattacharva & Kosowskv 2007l: IHemandez-Monteagudo et al 20061 : 
DePeo, Spergel & Trac 2005; iMaturi et al. 20071 : 



RoncarelU et al. 2007h . But one of the major sources of un- 



certainty in modeling the kSZ effect is the gas fraction and its 
evolution. So understanding both the thermal and kinematic SZ 
signals requires detailed understanding of feedback mechanisms 
in galaxy clusters and groups. The mechanisms and effects of 
feedback are also a long-standing question in astrophysics, with 
particular bearing on the process of galaxy formation. 

To this end, we have analyzed a sample of ten galaxy groups 
at z = 1 from numerical cosmological simulations of gas and 
dark matter which have been extended to include a self-consistent 
model for the evolution of massive black holes and their baryon 
feedback. At redshift z > 1, the quasar mode of black hole ac- 
cretion is expected to be the dominant feedback mechanism, com- 
pared to the radio-lou d accretion mode which becomes impor- 
tant at lower redshifts jSiiacki et al. 20071) . The size of our simu- 
lations prevents studying feedback in galaxy clusters, but rather re- 
stricts us to less massive galaxy groups. But as already mentioned, 
galaxy groups with shallow potential wells provide the best place 
to study non-gravitational heating and its implications for the prop- 
erties of hot gas. High-redshift galaxy groups are also a major con- 
tributor to the thermal SZ power spectrum, which peaks around 
z fa 1, when galaxy groups a re more numerous than massive clus- 
ters ( iKomatsu & Seliak 20021) . 

Following this introduction. Section II describes our simula- 
tion and its implementation of quasar feedback. In Section III we 
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Box size 


Np 
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■^cnd 








H-^Mq 




kpc 




D4 


33.75 


2 X 216=* 


2.75 X 10* 


4.24 X 10'^ 


6.25 


0.00 


D6 


33.75 


2 X 486=* 


2.75 X 10^ 


4.24 X 10*' 


2.73 


1.00 



Table 1. Numerical parameters of cosmological simulations (D4 & D6). 



study the effect of numerical resolution on our results; in Section IV 
we describe our results and compare them with a simulation that do 
not include quasar feedback. Finally, in Section V we summarize 
our results and discuss directions for future work, including mo- 
tivations and prospects for studying more massive galaxy clusters 
and more realistic feedback modeling for quasars and AGN. 



2 SIMULATION 

The cosmological simulati ons use d in this study are described 
in detail in tPi Matteo et al. 2007h . They use an LCDM cosmo- 
logical mode l with parameters c onsistent with the WMAP first- 
year results dSpergel et al. 20031) : fio ~ 0.3, JIa = 0-7, pri- 
mordial power spectral index n = 1, Hubble parameter h — 
0.7 with Ho = lOO/i kms~^ Mpc~^, and matter power spec- 
trum normalization erg = 0.9. A Gaussian random initial con- 
dition for this cosmology is evolved from high redshifts to the 
current epoch using a modified version of the parallel TreePM- 
SPH code GADGET2 jSpringel 20o3) . which manifestly con- 
serve entropy and energy. Gas dynamics is implemented with 
the L agrangian smooth ed-particle hydrodynamics (SPH) tech- 
nique JMonaghan 1992h . Radiative cooling and heating processes 
are co mputed with a spatially uniform ph otoionizing UV back- 
ground ( iKatz. Weinberg & Hemguist 19961) . For modeling star for- 
mation and its associated supernova feedback the code uses a sub- 
res olution multiphase model for the interstellar medium developed 
by jSpringel & Hemguist 20030 . In this model, a thermal instabil- 
ity is assumed to operate above a critical density threshold pth, 
producing a two phase medium consisting of cold clouds embed- 
ded in a tenuous gas at pressure equilibrium. Stars form from the 
cold clouds, and short-lived stars supply an energy of lO^'^ ergs 
to the surrounding gas as supernovae. This energy heats the dif- 
fuse phase of the ISM and evaporates cold clouds, thereby estab- 
lishing a self-regulation cycle for star formation. The pth is deter- 
mined self-consistently in the model by requiring that the equation 
of state (EOS) is continuous at the onset of star formation. The 
cloud evaporation process and the cooling function of the gas then 
determine the temperatures and the mass fractions of the two hot 
and cold phases of the ISM, such that the EOS of the model can 
be directly computed as a function of density. The latter is encap- 
sulating the self-regulated nature of star formation owing to su- 
pem ovae feedback in a simple m odel for a multiphase ISM. As in 
the ( Springel & Hernguist 2 0031 model we have included a model 
for supernova-driven galactic winds with an initial wind speed of 
V ~ 480kms~\ 

A unique aspect of the simulations is their inclusion of su- 
permassive bla ck holes and the resu lting energy feedback from 
mass accretion Matteo et al. 20071) . Black holes are represented 
as collisionless "sink" particles which grows from a seed black 
hole through accretion of mass from its immediately surrounding 
gas or through merger with another black hole. Seed black holes 
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Groups -R200m -RsOOm M200m -M'sOOm 

Mpc/h Mpc/h IQI^Mq/Ii IQl^Mg/h 






0.80 


0.56 


4.71 


3.08 


1 


0.77 


0.57 


4.40 


3.10 


2 


0.75 


0.45 


2.97 


1.57 


3 


0.68 


0.46 


2.14 


1.64 


4 


0.65 


0.41 


1.89 


1.21 


5 


0.63 


0.36 


1.780 


0.82 


6 


0.63 


0.37 


1.783 


0.84 


7 


0.60 


0.36 


1.47 


0.80 


8 


0.57 


0.34 


1.23 


0.67 


9 


0.53 


0.36 


1.13 


0.76 



Table 2. Properties of galaxy groups in the simulations at 2 = 1 

of mass M — l(f h^^ Mq are placed into the centers of ha- 
los whenever they reach a mass threshold of 10"'/i~^Mq. The 
subsequent gas accretion rate onto the black hole is estimated 
using the Bondi-Hoyle-Lyttleton parametriz ation jSondi 19521 ; 
iBondi & Hovle 1944l : Hovle & Lvttleton 19391) . We assume a fixed 
value 77 = 0.1 for the radiative efficiency rj = Lr/{MBHC^), 
where Lr is the radiated luminosity and Mbh is the mass accre- 
tion rate. This efficiency value is the mean value of a radiatively 
efficient accretion disk onto a Schwarzschild black hole Q). We 
further assume that a fraction e/ of i,, couples to the surround- 
ing gas in the form of feedback energy Ef deposited isotropi- 
cally, i.e. Ef = e/I/r- A fixed value of e/ = 0.05 is adopted 
here to fit current data on the normalization of the A/bh — cf 
relation between black hole mass and stellar velocity dispersion 
loi Matteo. Spri ngel & Hemquist 2005h . 

We use three different simulation runs, each of box size 
33.75 Mpc//i. The box size is a compromise between the re- 
quirements of sufficient spatial resolution to resolve physical pro- 
cesses in high-density regions surrounding black holes and suf- 
ficient volume to allow formation of halos with galaxy group 
masses. We study halos at z = 1: below this redshift, the funda- 
mental modes in the cosmological box become nonlinear and the 
simulations become u meliable on scales of their largest objects 
jPi Matteo et al 2003h . We name the runs D4 (with and without 
black holes) and D 6 (include black holes) foll owing the naming 
scheme adopted in jSpringel & Hemgiiist 2003h . Runs D4 and D6 
include black hole accretion along with cooling, star formation and 
supernova feedback, while the run-D4 (no black holes) leaves out 
black holes but includes all other physical processes. We use D4 
(no black holes) as a baseline comparison simulation to analyze the 
effects of quasar feedback on galaxy groups for the run D4. We also 
compare D4 and D6 to understand the issues of resolution and con- 
vergence. The numerical parameters of the runs, including particle 
number and mass resolution, are listed in Table[T] 

Table |2] lists the radius and mass of the galaxy groups formed 
in these simulations (the bulk group properties are essentially the 
same for all three simulations). Masses are defined as the amount 
of matter contained within a spherical region of overdensity 200 
(M2oo m) or 500 (Msppm) tim es the mean density of the universe at 
Matteo et al 20031) . Figure[T] shows gas density and star 
density for the most massive halo(A/2oom = 4.7 x lO'^/i-'M©) 
in the simulation. The left panel shows the map for each of the 
properties when black hole feedback is included while the right 
panel gives the map with no quasar feedback. Note the gas density 
maps are color coded by temperature- the brightness shows the den- 
sity and the color represents the temperature. It is evident that the 



gas is hotter when the feedback is included compared to when not 
included. Also the distribution of stars has changed significantly 
when quasar feedback is included. 



3 EFFECTS OF NUMERICAL RESOLUTION 

To study the effect of quasar feedback, we need to resolve quasars 
and surrounding gas at kiloparsec scales while simultaneously fol- 
lowing the formation and evolution of galaxy groups at megaparsec 
scales. Given this huge dynamic range, it is worthwhile to check 
how numerical resolution affects our results. We have run two 
simulations, namely "D4" and "D6," with the same cosmological 
parameters, initial conditions and simulation volume. The lower- 
resolution D4 run uses 2 x 216'^ total particles, while D6 uses 

2 X 486'^ particles. The corresponding mass resolution of the gas is 
4.24 X IO'^Mq /h and 4.24 x IQ^Mq /h. Their spatial resolution is 
characterized by gravitational softening lengths of 6.25 kpc/Zi and 
2.73 kpc//i respectively. 

We have studied the difference in the star and gas distributions 
at redshift z = 1, with comparisons displayed in Fig.|2] These plots 
show the average differential profile in the simulations. 

On average, both star and gas distributions agree within 10% 
for the D4 and D6 runs for R > 0.1i?200m - Beyond 7? = i?200m, 
statistical fluctuations causes star distributions to vary. Note that 
most of the star formation occurs in the inner region of the halo, so 
these statistical variations in the outer parts do not affect any of the 
conclusions about star fraction. 

The temperature profile shows roughly 10 — 15% difference 
between the simulations D4 and D6 in the inner region of the 
cluster, dropping to 5% for R > 0.2f?200m- The pressure profile 
shows relatively more robustness to numerical resolution: a 10%- 
15% difference in the inner region drops to only 5% to 3% for 
R > 0.2i?200m- Numerical resolution should thus have a minimal 
effect on the thermal SZ flux, since most of the signal comes from 
outside the core. Finally, the entropy profile shows a difference of 
20% in the inner region and a 5% diff erence for R > 0.3i?200m ■ 

As already show n in jSpringel & Hemquist 20031: 
iHemquist & Springel 20030 . using a large number of cosmo- 
logical simulations, simulation including star formation and 
cooling converge reasonably well in the resolution range between 
D4 and D6. 

Given this rough quantification of the effect of increased res- 
olution, we proceed to analyze the lower-resolution D4 simulation 
in the rest of the paper and compare with the same resolution run 
without black holes, noting where errors due to numerical resolu- 
tion limits might be a significant fraction of the effects being dis- 
cussed. In the following sections, we study the differential and cu- 
mulative profile for each physical quantity with two lower resolu- 
tion D4 runs both with and without including black holes. For each 
physical quantity we also calculate the difference in the profiles for 
each halo between the two runs with and without black holes and 
then show the mean of the difference. Also we find there are atleast 

3 mergers namely 2nd, 5th and 8th most massive halos(in the group 
of 10 halos) we considered here. While studying the average pro- 
files, we have excluded these halos from the averaging process so 
that the profiles do not get biased. However we have reported prop- 
erties of all the halos when they are studied as a function of mass. 
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no BH feedback 




20 kpc 



stars 



w BH feedback stars 




Figure 1. The gas distribution (top) and stai" distribution (bottom), both with quasar feedback (left column) and without (right column), for a halo of mass 
M = 4.6 X 10^^ A'/q at z = 1. The gas density maps are color coded by temperature (brightness shows density and color represents temperature). Note the 
qualitative difference in the distribution of stars between the two simulations. 



4 RESULTS 

4.1 Thermodynamics of the Intracluster Medium 

In this Section we study the impact of quasar feedback on thermo- 
dynamics of the ICM, namely on the three quantities pressure, tem- 
perature and entropy. Figure [3]gives the average temperature profile 
with scatter around the mean. In the inner region (i? < 0.2i?200m) 
of the halo, the temperature is enhanced by about 15-20% and by 
5-10% in the region 0.2_R200m < R < 0.57?2OOm. This is physi- 
cally reasonable as quasar feedback is coupling part of its radiated 
thermal energy to the surrounding ICM. We do not see any change 
in temperature due to quasar feedback at radii outside the halo core. 



For comparison we also show the average mean profile from the D6 
run. 

Note however, that the temperature profile inside the halo core 
becomes steeper when the feedback is included, whereas the obser- 
vations at low redshift shows a rather flat profile inside the core. 
This disagreement might be either due to the inability of the feed- 
back mechanism to explain the observed temperature profile and 
an improved model is needed or that one needs to include other 
sources of feedback in the simulations. Observations of group size 
halos at high redshift will be needed in order to understand whether 
the temperature profile indeed gets steeper at higher redshift or a 
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Figure 2. The effect of numerical resolution on vaiious quantities as functions of radius: baryons, both gas and stars (top left ), temperature (top right), pressure 
(bottom left) and entropy (bottom right). In each panel, dotted lines represents higher resolution (D6) and solid lines represent lower resolution(D4). 



better feedback mechanism is required to explain the flatness of the 
temperature profile. 

The temperature of the system agrees fa irly well with previ- 
ous studies made using halos of similar mass jBorgani et al. 2()04l : 
iKhalatvan et al. 20071 ; iFinoguenov et al. 200 ih . For eg., a halo of 
mass 4.7 X M^hT^ is expected to have a temperature of 
around 1 keV at z=0. We find a temperature of 1 .5 keV for a similar 
system at z=l. If a virial scaling relation is assumed this translates 
to a temperature of about 1 keV at z=0 which is consistent with 
previous studies. 

The corresponding average pressure profile is shown in Fig.|4] 
We find that the pressure decreases for R < 0.3 Mpc/h, beyond 
which quasar feedback clearly leads to a pressure enhancement of 



15% to 20% out to radius of R2oom- The entropy profile is shown 
in Fig.[5] The excess entropy near the core region is 50% larger than 
the no feedback case. The observational fin ding for the entropy pro- 
file for small groups jPonman et al. 20(331) agrees fairly well with 
the current study when virial scaling is assumed to translate the en- 
tropy profile at z=l in the current study to z=0. The scatter around 
the mean profile for each of these quantities is large, so we need a 
larger sample size to confirm these systematic deviations. The en- 
tropy and pressure profile indicates that the quasar feedback has 
driven the gas out from the inner region and redistributed in the 
outer region. The lower panels of the figures show the fractional 
difference for each quantity. As shown, in the inner region the dif- 
ference in the profiles is significant; far in excess of the numerical 
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Figure 3. The mean differential (left) and cumulative (right) temperature profile of gas averaged over seven halos. For each top panel, solid lines represent the 
mean and scatter around the mean profile for simulation D4 including quasar feedback, while the dotted lines represents the same quantities for simulation D4 
with no quasar feedback. Also shown is the mean profile from the D6 run (blue dashed line. The lower panels show the mean fractional change between the 
halos in the two runs. The blue dashed line shows the mean and the scatter in the difference in the profiles between D4 and D6 ( resolution effect) while the 
soUd red Une shows similar difference between the D4 runs ( the effect of including the black holes) 




Figure 4. Same as in Fig. [5] except for pressure. 



resolution error. Similar differences can be seen in the outside re- 
gion where the numerical resolution error is few percent. 



4.2 Baryon Fraction of the Intracluster Medium 

A particularly important issue for interpreting future Sunyaev- 
Zeldovich measurements is the gas fraction in a given halo. Here 
we consider the effect of quasar feedback on both baryonic com- 
ponents, stars and hot gas. The ten most massive objects formed 
in the simulation have masses ranging from 1 to 5 X 10^^ M©//;.. 
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Each object is binned in spherical shells, and the mass fractions of 
stars, gas and dark matter within each shell are normalized to the 
primordial baryon fraction Q,b/0.m,- Figure[6]shows the average dif- 
ferential (left) and cumulative (right) distribution of gas and stars. 
Note the difference in star formation between the simulations with 
and without quasar feedback is on average 20% to 40% out to ra- 
dius 7? = 0.67i200m- It is evident that quasar feedback substantially 
suppresses star formation at all radii; the cumulative star distribu- 
tion is 30% lower when feedback is included. The feedback mech- 
anism provides enough pressure support that a significant amount 
of gas fails to collapse and form stars. Comparing differential and 
cumulative profiles, it is evident that most of the star formation is 
suppressed in the interior region of the halo. 

Quasar feedback has an equally significant effect on the gas 
distribution. As shown in the top panel of Fig. |6l hot gas is being 
driven out from the internal region of the halo {R < Rsoom) to- 
wards the outer region. The gas density is lowered by 20-30% in the 
core; to compensate for this depletion, gas density is 10% higher at 
R > 0.3ii2oom compared to the no-feedback case. As is evident 
from the cumulative gas distribution, the feedback is not powerful 
enough to drive the gas from gravitational well of the halo. Note 
that there is still a difference in total gas mass of around 4% within 
a radius of 2i?200m which compensates for the lower star formation 
in these halos. 

Figure|7]shows cumulative gas and star fractions as a function 
of halo mass, measured out to radii i?200m, -Rsoom, and i?2500m, 
and also between R^oom and i?2500m ■ Table 3 to 6 gives the frac- 
tions for individual halos at these radii and also the mean and scat- 
ter. On average, cumulative star fractions shows a 30% depletion 
at all radii < Rsoom in simulation with quasar feedback; Gas frac- 
tions shows only mild change at 7?200m and i?5oom, although at 
J?2500m the gas fraction is about 15% lower with quasar feedback. 
When halo cores are excluded (i.e. between -Rsoom and i?2500m), 
the gas fraction is enhanced by about 10% in simulation with feed- 
back. This again shows that gas is driven off from the inner region 
of the halos to outer region. The gas fraction < Rsoom displays a 



Mass 



Quasar Feedback No Quasar Feedback 





fgas 


fstar 


fgas 


fstar 


4.71 


0.81 


0.16 


0.76 


0.21 


4.40 


0.79 


0.13 


0.78 


0.16 


2.97 


0.81 


0.13 


0.78 


0.19 


2.14 


0.80 


0.12 


0.82 


0.12 


1.89 


0.79 


0.12 


0.76 


0.17 


1.78 


0.77 


0.15 


0.79 


0.21 


1.78 


0.78 


0.11 


0.74 


0.16 


1.47 


0.77 


0.14 


0.76 


0.18 


1.23 


0.65 


0.12 


0.75 


0.19 


1.13 


0.76 


0.14 


0.66 


0.31 


Mean 


0.77 


0.13 


0.76 


0.19 


Scatter 


0.21 


0.13 


0.20 


0.21 



Table 3. Cumulative fractions of gas and stars out to R200m , both with and 
without quasar feedback. 



slight trend with mass in both simulations, although the star frac- 
tion shows no such effect. 

Figure [S] shows the cumulative ratio of gas to 
stars. This quantity plays an important role for deter- 
mining the cosmic matter density dWhite et al. 19931; 
Evrard 19 9?!; I Allen. Schmidt & Fabian 20021; 



lEttori et al. iool ; 

[Allen et al. 20041) . Usually it is assumed that this rat io is fixed a t 
any radius with negligible redshift evolution (Ettori et al. 20061) . 
we find that this assumption does not hold for either of the simula- 
tions. Without quasar feedback, the gas mass to stellar mass ratio 
changes roughly from 2 to 5, a factor of 2.5, between 0.3i?200m 
and 7?200m; for the simulation including quasar feedback the 
corresponding change in the ratio is slightly larger, from 2 to 7.5, a 
factor of 3.5. The ratio rises more steeply for the simulation with 
feedback and continues increasing beyond R200m ■ 
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Figure 6. Same as in Fig. [3] except for gas density (top panels) and star density (lower panels). 



4.3 Thermal Sunyaev-Zeldovich Decrements 

The thermal Sunyaev-Zeldovich distortion from quasar 
feedback has been stu di ed previously (see, e.g., 
jChatteriee & Kosowskv 20071 : IScannapieco et al. 20081) and 
references therein). This effect has a systematic impact on galaxy 
group-sized halos. As discussed above, the inaccuracy in the 
pressure profile due to numerical resolution limitations is on 
the order of 10% for R < 0.1-R200m, so we exclude the halo 
core region when calculating SZ distortions. This does not sub- 
stantially affect any of our results since the major contribution 
to the SZ signal comes from the region outside the halo cores 
jKomatsu & Seliak 20021) . We calculate the mean Coinpton 
y-distortion, which we denote Y, by integrating the gas pressure 



along the line of sight for each halo out to a radius of /?200m 
and over the projected cross-section of the cluster in comoving 
coordinates. Figure |9] shows the Y versus mass for the halos 
considered here, both with and without quasar feedback, with the 
lower panel showing the fractional change in Y . The individual 
halo F-parameters are given in Table |7] On average, the Y 
parameter changes by 6% (excluding the mergers) due to quasar 
feedback in these galaxy groups. Note that Y parameter in the run 
with the feedback shows both increase and decrease compared to 
the no feedback run as a function of mass. 

We also give a power law fit to the F- mass relation of th e 
form Y/E{zf^'-^ = 10''(A/2oom/10"MQ)°' jSehgal et al. 2007h . 
where a and /3 are fitting parameters and E{z) — + z)^ + 
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Mass (1013 h-iMj Mass (10>3 h-'Mj 



Figure 7. Cumulative gas and star fractions for the 10 most massive groups at z = 1 measured within a radius R = R200m (top left), R = RsoOm (top 
light), and R = R2500m (lower right), and between _R = -Rsoom and R = R2500m (lower left). For each panel, squares represents the star fraction and 
triangles the gas fraction. Solid lines correspond to the simulation including quasar feedback and dotted lines represent the no-feedback case. 



JIa)"'^ is the redshift evolution of the Hubble parameter. Although 
the scatter is large, the power-law fits in both simulations (given in 
Table [8] are close, and the val ues are consisten t with other studies 
with larger numbers of halos dSehgal et al. 20071) . 

As shown in jKomatsu & Seliak 2002h . the SZ power spec- 
trum receives a dominant contribution from high redshift halos; es- 
pecially for / > 3000, the contribution to Ci comes mostly from 
z > 1. The halo mass range considered here provides significant 
contribution to the Ci for / > 5000 and non-negligible contribution 
for / = 3000 to 5000. Since C; oc Y^, we expect that quasar feed- 



back will lead to a systematic increase in C; on the order of 10% 
between I = 5000 and 10000. 

Note that the difference in Y between the feedback and no- 
feedback cases does not tend to decrease with mass (Fig. |9}, al- 
though the scatter is too large to claim any statistical significance of 
this behavior. It is imperative to simulate bigger volumes to quan- 
tify the effect of quasar feedback on the V-mass relation for galaxy 
clusters, and the corresponding systematic differences in cluster 
mass estimates. We also emphasize that effect of quasar feedback 
generally increases with redshift, so our results at z = 1 give con- 
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Mass Quasar Feedback No Quasar Feedback 





fgas 


fstar 


fgas 


fstar 


4.71 


0.70 


0.21 


0.73 


0.31 


4.40 


0.77 


0.15 


0.75 


0.22 


2.97 


0.73 


0.19 


0.69 


0.28 


2.14 


0.69 


0.17 


0.72 


0.19 


1.89 


0.60 


0.18 


0.68 


0.30 


1.78 


0.54 


0.25 


0.72 


0.27 


1.78 


0.75 


0.14 


0.71 


0.19 


1.47 


0.64 


0.21 


0.71 


0.29 


1.23 


0.44 


0.17 


0.66 


0.29 


1.13 


0.71 


0.19 


0.59 


0.43 


Mean 


0.66 


0.18 


0.70 


0.28 


Scatter 


0.31 


0.18 


0.20 


0.25 



Table 4. Same as in Table[3] for -R500; 



Mass Quasar Feedback No Quasar Feedback 



10^=^ Mq/i-1 


fgas 


fstar 


fgas 


fstar 


4.71 


0.49 


0.34 


0.65 


0.45 


4.40 


0.61 


0.22 


0.66 


0.31 


2.97 


0.63 


0.27 


0.63 


0.43 


2.14 


0.65 


0.20 


0.69 


0.22 


1.89 


0.48 


0.22 


0.61 


0.40 


1.78 


0.49 


0.32 


0.63 


0.33 


1.78 


0.70 


0.17 


0.70 


0.28 


1.47 


0.55 


0.27 


0.70 


0.35 


1.23 


0.34 


0.22 


0.60 


0.35 


1.13 


0.65 


0.22 


0.56 


0.48 


Mean 


0.56 


0.24 


0.64 


0.36 


Scatter 


0.11 


0.054 


0.046 


0.081 



Table 5. Same as in Table[3] for -R2500m- 

servative estimates of the quasar feedback impact on the SZ signal 
at earlier times. 



5 DISCUSSION 

We have studied the effect of quasar feedback on baryon frac- 
tions and on thermodynamics of intracluster medium of interme- 
diate mass halos corresponding to galaxy groups. Our analysis uses 
high-resolution N-body plus hydro cosmological simulations in a 
box with side length 33.75 Mpc/ft. One simulation is conventional, 
while another incorporates black hole growth, accretion, and en- 
ergy ejection assuming simple astrophysics consistent with obser- 
vations; both simulations use the same initial conditions so indi- 
vidual large haloes can be compared. From the ten most massive 
haloes, with masses ranging between 1 and 5 x IO^^Mq/Zi, we 
draw the following conclusions: 

1. Compared to the no-feedback case, star formation is sup- 
pressed by 30-40% in the iimer regions of the halos because of the 
additional pressure support provided by quasar feedback. 

2. Quasar feedback redistributes hot gas, driving it from the 
inner region towards the outer part of the halos. As a result, gas 
density is 20% less in the inner part and 10% to 15% greater in the 
outer region when compared to the simulation without feedback. 
However, the gas fraction in the two simulation differs by only 5% 




Figure 8. The average cumulative fraction of the ratio of gas mass and 
stellar mass for the ten halos. Solid lines represent the mean and scatter for 
the sample including quasar feedback, while dotted lines represent the same 
for the no-feedback case. 



T 1 r 




5 

Mass (10'3 h-'Mj 

Figure 9. Sunyaev-Zeldovich V-distortion versus halo mass for 10 halos, 
for mass and gas within lMpc//i of the halo center. Squares represent 
values from simulation D4 including quasar feedback; triangles represent 
values from simulation D4 without feedback. Lines are the best-fit power 
law to the y-mass relation including quasar feedback (solid) and without 
quasar feedback (dotted). The lower panel shows the fractional change in 
Y between the two simulations. 
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Mass Quasar Feedback No Quasar Feedback 





fgas 


fstar 


fgas 


fstar 


4.71 


0.80 


0.11 


0.75 


0.31 


4.40 


0.83 


0.10 


0.76 


0.16 


2.97 


0.78 


0.12 


0.70 


0.29 


2.14 


0.71 


0.13 


0.73 


0.17 


1.89 


0.68 


0.09 


0.68 


0.30 


1.78 


0.59 


0.18 


0.72 


0.27 


1.78 


0.82 


0.07 


0.70 


0.20 


1.47 


0.73 


0.09 


0.73 


0.16 


1.23 


0.51 


0.09 


0.66 


0.30 


1.13 


0.79 


0.08 


0.61 


0.45 


Mean 


0.72 


0.11 


0.70 


0.26 


Scatter 


0.10 


0.032 


0.04 


0.09 



Table 6. Cumulative fractions of gas and stars between radii Rsoom and 
^J2500m, both with and without quasar feedback. 



M200m 


Y (feedback) 


Y (no feedback) 






10-'^ Mpc2 


lO-^Mpc^ 


% 


4.71 


1.91 


1.88 


0.02 


4.40 


1.43 


2.06 


-0.44 


2.97 


0.71 


0.67 


0.04 


2.14 


0.54 


0.48 


0.12 


1.89 


0.37 


0.37 


0.01 


1.78 


0.26 


0.25 


0.03 


1.78 


0.18 


0.24 


-0.33 


1.47 


0.21 


0.24 


-0.11 


1.23 


0.19 


0.21 


-0.08 


1.13 


0.18 


0.16 


0.11 



Table 7. The relation between SZ Y-distortion and cluster mass for galaxy 
groups with and without quasar feedback. 

to 10%, and gas fractions tends to increase mildly with increasing 
halo mass. 

3. The ratio of gas mass to stellar mass increases by a fac- 
tor of 3.5 in the simulation including quasar feedback and a fac- 
tor of 2.5 in the simulation without quasar feedback in the region 
0.2ii200m < R < 0.5-R200m- This contradicts the common as- 
sumption that this ratio is constant at all radii. 

4. Both temperature and entropy increase by 30% to 50% in 
the halo core region because of the additional thermal energy radi- 
ated by quasars. 

5. Pressure decreases by 30% in the inner region and increases 
by 15% to 20% at radii larger than 0.4 i?200m due to the increased 
gas density in this region. This leads to a change of about 6% in the 
mean Sunyaev-Zeldovich V-distortion. The resulting SZ angular 
power spectrum will be larger by around 10% for I > 5000. We 
find little dependence of the SZ enhancement with halo mass. 



a 


/3 


with feedback 1.78 ±0.06 


-5.55 ±0.17 


no feedback 1.79 ±0.05 


-5.47 ±0.13 



Table 8. Power law fits to the SZ y-mass relation for galaxy groups with 
and without quasar feedback, as displayed in Fig.|9] 



The effects of quasar feedback on the intracluster medium will 
be most evident in the group-sized haloes considered here, with 
their relatively shallow gravitational potential wells. Observation- 
ally, the most interesting haloes are larger in mass by a factor of 
ten, galaxy clusters: these are the haloes which are most readily de- 
tected via their SZ, X-ray, or optical signals. The gas fractions do 
not show any particular trend with increasing halo mass, and star 
fractions increase very weakly with mass over the halo mass range 
studied here. So it is reasonable to expect that the results of this 
work will hold for cluster-sized halos as well. Nevertheless, given 
the substantial impact of quasar feedback on various properties of 
the intracluster medium which the current study suggests, it is im- 
perative to study cluster-sized halos as well. This requires larger- 
volume simulations, as the number density of clusters decreases 
with cluster mass. To this end, we are currently running a simula- 
tion of box size 50 Mpc/^; results will be reported elsewhere. 

This is the first attempt to study the impact of quasar feed- 
back on the baryon fraction and thermodynamics of the intraclus- 
ter medium in a cosmological hydrodynamic simulation. Both the 
gas and star fractions in our simulation are consistent with cur - 
rent observational limits l l Allen et aL 20041: lEttori & Fabian 19991) . 
Note that we have studied only quasar feedback at redshifts greater 
than unity. However, active galactic nuclei also inject energy into 
the ICM via a "radio mode" which is be lieved to be the dominant 
feedback mechani sm at lower redshift jSiiacki & SpringeT 20061: 
ISiiacki et al. 2007|) . Thus our results should be treated as a con- 
servative estimate of the total impact of AGN feedback for galaxy 
groups at low redshifts. 

Gas pressure in cosmological halos, particularly those with 
masses ranging from galaxy groups to galaxy clusters, determines 
the important thermal Sunyaev-Zeldovich signal which will soon 
be measured with high precision. The gas fraction is important for 
connecting kinematic SZ signals of cluster gas momentum with 
theoretical predictions about cluster velocity or total momentum. 
This paper takes the first step towards quantifying the impact of 
quasars on these quantities, which turns out to be significant but not 
dominating. Much work remains to be done, both through larger 
simulations which contain many galaxy-cluster-sized haloes, and 
in enhancing the realism of the quasar feedback models. We hope 
the results here plus the exciting observational prospects in the near 
future will open the door to further advances in this area. 
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